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Abstract 

We calculate accurate eigenvalues of a bounded oscillator by means 
of the Riccati-Pade method that is based on a rational approximation 
to a regularized logarithmic derivative of the wavefunction. Sequences 
of roots of Hankel determinants approach the model eigenvalues from 
below with remarkable convergence rate. 

1 Introduction 

Some time ago Fernandez et al [1-11] developed a method for the accurate 
calculation of eigenfunctions and eigenvalues for bound states and resonances 
of separable quantum-mechanical models. This approach is based on Pade 
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approximants built from the Taylor expansion of a regularized logarithmic 
derivative of the eigenfunction. If such rational approximations are forced 
to yield an additional coefficient of the Taylor expansion, then the physical 
eigenvalue is given by a convergent sequence of roots of Hankel determinants 
constructed from the coefficients of that series. One merit of this approach, 
called Riccati-Pade method (RPM), is the great convergence rate of such 
sequences in most cases, and that the same equation applies to bound states 
and resonances [1-11]. Besides, in some cases the sequences yield upper and 
lower bounds to the eigenvalues [2]. 

The approach just mentioned was mainly applied to quantum-mechanical 
models with asymptotic boundary conditions at infinity [1-11], but it was 
suggested that it may be suitable for bounded problems as well [8]. In this pa- 
per we try the RPM on the bounded anharmonic oscillator recently discussed 
by Ciftci and Ateser [12] that may be suitable for the study of quarkonium 
physics. 

In Sec. [2] we outline the main features of the RPM. In Sec. [3] we apply the 
RPM to the bounded oscillator proposed by Ciftci and Ateser [12]. In Sec. H] 
we briefly address the problem of spurious roots. In Sec. [5] we show that the 
RPM applies to the central-field version of the bounded oscillator. Finally, 
in Sec. |6] we draw some conclusions. 

2 Method 

In this section we outline the main features of the Riccati-Pade method 
(RPM). For concreteness we consider the dimensionless Schrodinger equation 
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with an even-parity potential V(x) 

ip"(x) + [E — V(x)]ip(x) = (1) 

with yet unspecified boundary conditions. We assume that the Taylor ex- 
pansion 

oo 

V(x) = £ V j3 * (2) 

3=0 

converges in a neighbourhood of x = 0. 

The modified logarithmic derivative f(x) = s/x — ip'(x)/ip(x) satisfies the 
Riccati equation 



f'{x) + -f(x) - f(x) 2 - + V(x) -E = (3) 

X X A 



where the fourth term disappears when s = or s = 1 that correspond to 
even or odd states i/j(x), respectively. The Taylor expansion 
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/(•'■) •''X/r'' 2 ' (4) 

3=0 

also converges in a neighbourhood of x = and the coefficients fj depend on 
the eigenvalue E. 

The Pade approximant 



^2j=0 &jX 



l N + d / N \ = ■'■ t'v" , (5) 



yields exactly the first 2N+d+l coefficients of the series (jl]). We look for the 
values of E that allow the Pade approximant to produce the next coefficient 
f2N+d+i exactly. Such energy values are given by the roots of the Hankel 
determinant Hfy = with matrix elements fd+i+j-i, i,j = 1,2, . . . , D, where 
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D = N + 1 is the matrix dimension, and d = 0, 1, . . . is the displacement 
[1-11]. 

Each Hankel determinant exhibits many roots but one easily identifies 
convergent sequences E\^\ D = 2, 3, . . . that approach the actual eigenvalues 
E n of the Schrodinger equation ([I]). If one such sequence is monotonously 
increasing (decreasing), then each element provides a lower (upper) bound 
to the eigenvalue. This property of the Hankel sequences was already proved 
for some problems [2]. 

Another interesting feature of the RPM is that one does not have to take 
into account the boundary conditions explicitly. The method selects the 
appropriate values of E according to the singularities of the potential as we 
shall see below. 



3 The model 

Ciftci and Ateser [12] studied the Schrodinger equation ([1]) with the potential 

^2 2 

V(x) = — . (6) 

(1 -x 2 /R 2 ) 2 

where a,R > 0, —R < x < R, and the boundary conditions if)(±R) = 
which are obviously determined by the poles of V(x). We begin our analysis 
by noticing that the change of variables x = Lq leads to 

$(g) = (7) 

where $(g) = ip(Lq). Therefore, it is clear that one can choose either a or R 
equal to unity without loss of generality as follows from the obvious equalities 
E(a, R) = R~ 2 E(aR 2 , 1) = aE(l, aR). From a purely numerical point of view 
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it is convenient to calculate R 2 E(1,R) for small R and E(1,R) for large R 
because \im R _ R 2 E n (l, R) — (n + l) 2 7r 2 /4 and lim^oo E n (l, R) — 2n + 1, 
where n = 0, 1, . . . (from now on it should be assumed that a = 1 unless 
otherwise stated). 

We first show that the RPM yields the exact harmonic oscillator eigen- 
values when R —> oo. For example, the determinant of smallest dimension 
for s = and d = is 

h0 _ E 6 R 4 - 27E 4 R 4 + 32AE 3 R 2 + 51E 2 R 4 - Q75E 2 - 32AER 2 - 25i? 4 - 81 
2 ~ 4725i? 4 

(8) 

so that 

lim (g »-25)(g-l)' 

i?-oo 2 4725 v ; 

which clearly shows that two roots of the Hankel determinant give exactly the 
eigenvalues for the first two states of even parity of the harmonic oscillator. 

On the other hand, present implementation of the RPM is unable to 
yield exact eigenvalues for R < oo, and the convergence rate of the sequence 
of roots of the Hankel determinants decreases as R decreases. However, the 
RPM provides remarkably accurate results for all R as shown in what follows. 

Table [T] shows two sequences of roots of the Hankel determinants that 
converge towards the ground-state eigenvalue of the bounded oscillator given 
as R 2 Eq(\, R), for R = 0.1. The first ten digits of present result agree with 
those reported by Ciftci and Ateser [12]. The velocity of convergence of the 
RPM sequences is as remarkable as in previous applications of the method 
to problems with asymptotic boundary conditions at infinity [1-11]. 

We have decided to keep 20 digits in all our results although in some cases 
not all of them are exact. In this way one can clearly estimate the accuracy 



of an entry by its agreement with the neighbours in the same column and 
row, making it easier to appreciate the appearance of new stable digits as D 
or d increases. 

Table [2] shows the ground-state energy for some values of R and the 
smallest determinant dimension at which all the reported digits are stable 
(d = in all cases). 

The convergence rate of the RPM sequence decreases as the quantum 
number increases. The reason is that N should increase in order to take into 
account the increasing oscillation of the excited states. Table [3] shows two 
RPM sequences converging towards the eigenvalue E2 for R — 1. The first 
sixt digits of our result agree with those reported by Ciftci and Ateser [12]. 

4 Spurious roots 

As said above, the Hankel determinants exhibit many roots and one has to 
select those corresponding to the chosen state and model. Equation (Q shows 
that there is a double root corresponding to the ground state E — 1, and 
a single root corresponding to the next higher even state E2 = 5. Multiple 
roots in an exactly solvable problem become multiple sequences converging to 
the eigenvalue of a nontrivial problem [1-11]. However, one easily identifies 
the sequence with the best convergence rate looking for the root of H^ +1 
closest to the previously chosen root of Hfy. As D increases, more roots 
appear in the neighbourhood of the actual eigenvalue which is an additional 
indication of sound convergence. 

The change of variables x —> ix transforms the Hamiltonian operator 
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with potential V(x) = wz^rnpyi into minus the Hamiltonian operator with 
potential V 2 (x) = ( 1+ Jy fi 2)2 ■ This transformation accounts for the spurious 
roots E — — 1 and E = — 5 in equation (jUJ) (for i? — > oo) as discussed 
earlier in a more general way [8]. When R < oo the potential V^Oe) does not 
support bound states (the RPM does no longer select Dirichlet boundary 
conditions at ±.R, and — oo < x < oo). In this case one expects resonances 
with complex values of E, and outgoing or ingoing waves as x — > ±oo [6]. 
For example, when R = 1 there is a sequence of Hankel roots that converge 
towards E = —0.015565600439810503080 that appears to be minus the real 
part of a resonance of Vi(x). The complex part of this resonance seems to be 
quite small. We do not discuss this feature of the RPM any further because 
we are mainly interested in bounded oscillators. 

5 Central— field models 

Ciftci and Ateser [12] mention a possible application of a three-dimensional 
bounded oscillator to quarkonium physics. The application of the RPM to 
central-field models is straightforward [1,3,5,6]. The Riccati equation for 
the radial part of the Schrodinger equation becomes 

f(r) + -f(r) - /(r) 2 - + V(r) + -E = (10) 

where f(r) = s/r — ip'(r) /ip{r) and ip(0) = 0. In this case I = 0, 1, . . . is 
the angular-momentum quantum number and we choose s — I + 1 in order 
to remove the centrifugal term and make f(r) regular at origin [1,3,5,6]. 
When I = the results are identical to those of the odd states of the one- 
dimensional problem that satisfy the same boundary condition at origin. 
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Table H] shows two Hankel sequences converging towards E\ (or the energy of 
the ground state of the central-field model) when R — 1. The first six digits 
of our result agree with those reported by Ciftci and Ateser [12]. 

It is worth mentioning that the sequences shown in tables [U [31 and H] are 
monotonously increasing, and, therefore, they provide lower bounds to the 
corresponding eigenvalues of this particular model. 

6 Conclusions 

We have shown that the RPM gives accurate eigenvalues for some one- 
dimensional and central-field bounded oscillators. The approach automati- 
cally selects the physical energies according to the singularities of the potential- 
energy function. In addition to enlarging the field of application of the RPM 
we have verified that the perturbation variant of the AIM proposed by Ciftci 
and Ateser [12] yields reasonably accurate results for most physical purposes. 
However, the RPM appears to be much more accurate and straightforward. 
It does not require factorization of the asymptotic behaviour of the wavefunc- 
tion, and its Hankel sequences converge monotonously from below, providing 
lower bounds to the eigenvalues of this particular model. In this respect the 
RPM is a suitable complement to the variational methods that render upper 
bounds. 

Besides, the AIM has not yet proved suitable for the calculation of com- 
plex eigenvalues that the RPM yields easily and accurately [5,6,8,9,11]. It 
is worth mentioning that the RPM is not restricted to the Schrodinger equa- 
tion. We have recently applied a variant of the RPM, which we may call 
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Pade-Hankel method, to nonlinear two-point boundary value problems, ob- 
taining very accurate results for the unknown parameters in several models 
of physical interest [13]. 

References 

[1] F. M. Fernandez, Q. Ma, and R. H. Tipping, Phys. Rev. A 39 (1989) 
1605-1609. 

[2] F. M. Fernandez, Q. Ma, and R. H. Tipping, Phys. Rev. A 40 (1989) 
6149-6153. 

[3] F. M. Fernandez, Phys. Lett. A 166 (1992) 173-176. 

[4] F. M. Fernandez and R. Guardiola, J. Phys. A 26 (1993) 7169-7180. 

[5] F. M. Fernandez, Phys. Lett. A 203 (1995) 275-278. 

[6] F. M. Fernandez, J. Phys. A 28 (1995) 4043-4051. 

[7] F. M. Fernandez, J. Chem. Phys. 103 (1995) 6581-6585. 

[8] F. M. Fernandez, J. Phys. A 29 (1996) 3167-3177. 

[9] F. M. Fernandez, Phys. Rev. A 54 (1996) 1206-1209. 
[10] F. M. Fernandez, Chem. Phys. Lett 281 (1997) 337-342. 
[11] F. M. Fernandez, J. Phys. A 41 (2008) 065202 (5pp). 
[12] H. Ciftci and E. Ateser, Phys. Scr. 76 (2007) 517-520. 



9 



[13] P. Amore and F. M. Fernandez, Rational Approximation for Two-Point 
Boundary value problems, arXi v:0705. 3862 



10 



Table 1: Ground-state energy of the bounded oscillator times R 2 for R = 0.1 



D 


d = 


d = 1 


3 




2.3697606944397752864 


4 


2.4662533354307466207 


2.4674298730367798888 


5 


2.4674514194162024875 


2.4674515378204140996 


6 


2.4674515390300695979 


2.4674515390347723203 


7 


2.4674515390348029376 


2.4674515390348030263 


8 


2.4674515390348030267 


2.4674515390348030267 


9 


2.4674515390348030267 


2.4674515390348030267 


10 


2.4674515390348030267 


2.4674515390348030267 



Table 2: Ground-state energy of the bounded oscillator for some values of R 

R Eq D 

1 2.8848849919939971927 7 

10 1.0150378624976100592 3 

100 1.0001500037503748711 2 
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Table 3: Energy of the first even excited state for R = 1 



D 


d = 


d=l 


4 


24.086798714692429504 


24.361407377724659906 


5 


24.424509979572052554 


24.428687415859914431 


6 


24.429124533510868656 


24.429142207076630442 


7 


24.429143334496913081 


24.429143367079353985 


8 


24.429143368495260376 


24.429143368526373359 


9 


24.429143368527356872 


24.429143368527373946 


10 


24.429143368527374357 


24.429143368527374363 



Table 4: Energy of the first excited state of the bounded oscillator with R — 1 



D 


d = 


d= 1 


3 


11.134110459186094281 


11.158800510018056742 


4 


11.161646658411540522 


11.161732939936591584 


5 


11.161737758990207575 


11.161737854670925624 


6 


11.161737857894423316 


11.161737857940626897 


7 


11.161737857941670476 


11.161737857941681848 


8 


11.161737857941682032 


11.161737857941682034 


9 


11.161737857941682034 


11.161737857941682034 


10 


11.161737857941682034 


11.161737857941682034 
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